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Abstract. We present a brief (mainly bibliographical) report on recently performed calculations of terms 
of order 0{a'^n'j) and 0{a'^n^fmg) for hadronic Z and r decay rates. A few details about the analytical 
evaluation of the masters integrals appearing in the course of calculations are presented. 



PACS. 12.38.Bx Perturbative calculations - 12.38.-t Quantum chromodynamics 



1 Introduction 



Important physical observables like the cross-section of 
e~'"e~ annihilation into hadrons and the decay rate of the 
Z boson are related to (the absorptive parts of) the vector 
and axial- vector current correlators (for a detailed review 
see, e.g. ^). From the theoretical viewpoint the two- 
point correlators are ideally suited for evaluations in the 
framework of perturbative QCD (pQCD) 2^. Indeed, due 
to the simple kinematics (only one external momentum), 
even multiloop calculations can be analytically performed. 

In many important cases (with the Z- and r-decay 
rates as prominent examples) the external momentum is 
much larger than the masses of the relevant quarks. It 
is then justified to neglect these masses in a first ap- 
proximation which significantly simplifies all the calcu- 
lations. Within this massless approximation of QCD, the 
absorptive parts of vector and scalar correlators are ana- 
lytically known to '^T'E'B' . The residual quarks mass 
effects can be taken into account via an expansion in quark 
masses. This has been done for the quadratic and quartic 
terms to the same order jJIBlF!- 

During the past years, in particular through the anal- 
ysis of Z decays at LEP and of r decays, an enormous re- 
duction of the experimental uncertainty has been achieved. 
Inclusion of the 0{a^) corrections is mandatory already 
now. Quark mass effects must be included for Z-decays. 
The remaining theoretical uncertainty from uncalculated 
higher orders is at present comparable to the experimental 
one Thus, the full calculation of the next contributions, 
those of 0{a'^), to the two-point quark current correlators 
is an important next step in testing the Standard Model 
and crucial for precise determination of the QCD coupling 
constant. 



* To appear in the Proceedings of the International Euro- 
physics Conference on High Energy Physics, Aachen, Germany, 
17-23 July 2003. 



In this paper we discuss some selected theoretical as- 
pects of recent calculations of terms of order 0{a^nj) 



and 0{ag 



contributing to the two-point correlator 



of the vector (and axial) quark currents. More detailed 
discussion as well as phenomenological applications in the 
context of hadronic Z and r decays can be found in the 
original publications |11)II11II12| . 



2 Correlator and diagrams 

A well-known definition of a two-point quark current cor- 
relator reads 



I / dxe^^^(T[jX/f(x)(jX{f)nO)]) 



(1) 



with m? 



Y.f^) and = qr^^,{'^z)qj. The two (not 

necessarily different) quark fields with masses and mj 
are denoted by qi and qj respectively. 

The number of diagrams contributing to ^ grows fast 
with the order of perturbation theory. While only three di- 
agrams appear at 0{as), this number becomes 37 and 738 
at 0{a1) and 0(0^) respectively^. Finally, one arrives at 
19832 five-loop diagrams at C(a^). Even more important 
is the fact that the calculational complexity of a single 
diagram also grows tremendously with every additional 
loop. 

The calculational effort for a full evaluation of i?(s) in 
C'(af) is enormous and with present techniques exceeds 
the available computer resources by an order of magni- 
tude. It is for this reason that our calculations were lim- 
ited to a gauge invariant subset, namely the terms of order 
a^rij, where nj denotes the number of fermion fiavours. 

^ The specific numbers are cited as produced by the diagram 
generator QGRAF 13 for the case of the non-diagonal quark 
current (i ^ j in 10). 
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Some typical representatives of the corresponding set of 
diagrams are depicted in fig. ^ Note that the terms of 




Fig. 1. Some representative of four-loop diagrams contribut- 
ing to 10 in the 0{atn'f) order. 

order a^n^ are rather simple. They had been obtained 
earlier by summing the renormalon chains jl4j . 

3 Reduction to masters 

As is well-known 15 the calculation of absorptive part of 
the correlator ^ in the massless limit^ at the L + 1 loop 
level is reducible to the evaluation of some properly con- 
structed set of massless L-loop propagator-type diagrams 
(the corresponding Feynman integrals will be referred to 
as "p-integrals"). The completely automatized reduction 
procedure of [1] is based the method of Infrared Rear- 
rangement enforced by a special technique of dealing 
with infrared divergences ~ the i?*-operation |17j . 

Thus, to compute the five loop ©(a^) contribution in 
the absorptive part of ^ one should be able to evaluate 
generic four loop p-integrals. Unfortunately, this problem 
is far more complicated than the one at three loops. The 
latter was in principle done long ago by a manual con- 
sideration of all possible cases within the integration-by- 
parts technique (IBP) JF. Nevertheless, it took almost 
ten years before the corresponding algorithm was reliably 
implemented in FORM |19II2()| . A straightforward exten- 
sion of the same approach to four loops is barely possible 
at all. 

Fortunately, another method has been developed in 



According the IBP paradigm every integral is to be 
reduced to a sum of irreducible ("master") integrals, with 
coefficient functions being rational functions of the space- 
time dimension D. However, in contrast to the standard 
approach, where these coefficients are calculated by a re- 
cursive procedure, they are in the current approach ob- 
tained from an auxiliary integral representation 21 in 
the form of an expansion in 1/D. Calculating sufficiently 
many terms in this expansion, the original D-dependence 
can be reconstructed. 

The calculations were done in the following way. First, 
the set of irreducible integrals involved in the problem 

^ The statement also holds for the quadratic in quark masses 
contribution. 




Fig. 2. A generic massless two-loop diagram whose expansion 
in e was computed to order in |26| 

was constructed, using the criterion irreducibility of Feyn- 
man integrals Second, the coefficients multiplying 
these master integrals were calculated in the 1/D ^ 
expansion. This part was performed using the parallel ver- 
sion of FORM 25 running on an 8-alpha-processor-SMP- 
machine with disk space of 350 GB. The calculations in 
the massless limit took approximately 500 hours in total. 
(Approximately the same time went into the calculation 
of the ©(afn^m^) contribution.) Third, the exact answer 
was reconstructed from results of these expansions. Ex- 
tensive tests were performed. 



4 Master integrals and their evaluation 

The master integrals appearing in the course of our cal- 
culations are shown in Table ^ These can be separated 
in three group: simple (mOl, mil, ml2, ml3, ml4, m23, 
m24, m25, m31), semi-simple (m22, m26, m27, m21, m32, 
m33) and the difficult ones (m34, m35,m41 and m52). 

Simple integrals can be immediately performed in terms 
of r functions by a repeated application of the textbook 
one-loop integration formula. The members of semi-simple 
group prove to be easily reducible to a generic two-loop 
integral shown on fig. [21 The latter was computed up to 
pretty high order in the parameter e = 2 — D/2 in the 
eighties j26M27i ■ Luckily enough, remaining four difficult 
integrals were all analytically evaluated"^ long ago in the 
course of the computation of the five-loop /3-function in 
the (^''-theory ^51!^ . 



5 Summary 

The calculations performed in |10II11II12| demonstrate that 
the approach based on the 1 /D-expansion is suited to ob- 
tain genuine QCD results in five loop approximation. Fur- 
ther analysis shows that the other diagrams appearing in 
order can be solved in the same way, given sufficient 
computer resources. Work in this direction is in progress. 

This work was supported by Sondcrforschungsbereich- 
Transregio "Computational Particle Physics" (SFB-TR 9), 
by INTAS (grant 00-00313), by RFBR (grants 01-02-16171, 

^ We mean the pole and constant parts of the corresponding 
e expansions. 
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03-02-17177) and by the European Union under contract 
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Table 1. Master four- loop propagator- type integrals appear- 
ing in the course of calculations of the correlator (0 in orders 
Oiatn}) and 0{ain}ml). 



m41 



m34 





m35 



m52 




ni33 



m32 





m27 



m26 





m22 



m21 





m25 



m24 




ml4 



ml3 





m23 



mil 



m31 



rnl2 



<xxx> 

mOl 
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